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Abstract 

We consider the non-supersymmetric models of chaotic (driven by a quadratic potential) and hybrid 
inflation, taking into account the minimal possible radiative corrections to the inflationary potential. We 
show that two simple coupling functions /(ct) (with a parameter c-ji involved) between the inflaton field 
a and the Ricci scalar curvature ensure, for sub-Planckian values of the inflaton field, observationally 
acceptable values for the spectral index, n s , and sufficient reheating after inflation. In the case of 
chaotic inflation we consider two models with large c^'s resulting to n s ~ 0.955 or 0.967 and tensor- 
to-scalar ratio r ~ 0.2 or 0.003, respectively. In the case of hybrid inflation, the selected f(<j) assists 
us to obtain hilltop-type inflation. For values of the relevant mass parameter, m, less than 10 6 TeV 
and the observationally central value of n s , we find c-jz — (0.015 — 0.078) with the relevant coupling 
constants A = k and the symmetry breaking scale, M, confined in the ranges (2 ■ 10~ 7 — 0.001) and 
(1 - 16.8) • 10 17 GeV, respectively. 

Published in Phys. Lett. B 692, 287 (2010) 

1 Introduction 

Non-minimal inflation (non-MI) HI i.e. inflation constructed in the presence of a non-minimal cou- 
pling between the inflaton field and the Ricci scalar curvature, 1Z, has gained a fair amount of recent 
attention (2^4]]. In particular, it is shown that non-MI can be realized within the Standard Model (SM) 
- or minimal extensions of it - provided the inflaton couples strongly enough to 1Z. The role of 
inflaton can be played either by the Higgs doublet either by a SM singlet coupled to Higgs. Although 
quite compelling, non-MI within SM suffers from (i) several computational uncertainties regarding the 
impact of the quantum corrections in the presence of such a strong non-minimal gravitational coupling 
and (ii) the ambiguity about the hierarchy between the cutoff scale of the effective theory and the en- 
ergy scale of the inflationary plateau |6j|7). Be that as it may, it would be interesting to examine if 
appropriately selected non-minimal gravitational couplings can have beneficial consequences - as for 
the reconstruction of the cosmic expansion history [El - for other well-motivated and rather natural 
models of inflation (for a survey see, e.g., Ref. I9l). 

Two such models are undoubtedly Chaotic (CI) [ 10] and Hybrid Inflation (HI) ifTTTl . In this paper 
we focus on the non-supersymmetric version of these models. CI driven by a quadratic potential pro- 
vides the simplest realization of inflation without initial-value problem and with quite interesting pre- 
dictions for the (scalar) spectral index, n s , and the scalar-to-tensor ratio, r. However, trans-Planckian 
inflaton-field values are typically required to allow for a sufficiently long period of inflation. Thus non- 
renormalizable corrections from quantum gravity are expected to destroy the flatness of the potential, 
invalidating thereby CI. On the other hand, HI - although can be accommodated with sub-Planckian 
values for the inflaton - suffers from the problem of the enhanced n s which turns out to be, mostly, 
well above the prediction of the fitting fT2l of the five-year results from the Wilkinson Microwave 
Anisotropy Probe Satellite (WMAP5) plus baryon-acoustic-oscillations (BAO) and supernovae (SN) 
data - for an up-to-date analysis of the problem of initial conditions within HI, see Ref. Ifl3l . 
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Note, in passing, that the introduction of supersymmetry (SUSY) and its local extension - super- 
gravity (SUGRA) - can alleviate the shortcomings of both models - see Ref. lfl4l for several resolu- 
tions to the problem of CI and Ref. lfl5Hl9l for proposals related to the disadvantage of HI. However, 
we have to accept that there is no direct experimental confirmation of SUSY until now. On the con- 
trary, there is a strong observational evidence in favor of the inflationary paradigm. Consequently, it is 
worthwhile to build models of CI and HI consistently with the observations, even without the presence 
of SUSY - for similar recent attempts, see Ref. M201I21I1 . 

In this paper, we propose two types of non-minimal coupling functions f(a) between the inflaton 
and 1Z which support a resolution to the aforementioned problems of CI and HI. After the end of non- 
Mi, both /(<t)'s shrink to unit assuring thereby, a safe transition to the Einstein gravity in time. In the 
case of non-minimal CI (non-MCI), two models with clearly distinctive results are investigated. In the 
case of non-minimal HI (non-MHI), the inflationary trajectory is concave downwards and so, inflation 
turns out to be of hilltop-type 11221 . In both cases, the minimal possible radiative corrections |[23l 
to the inflationary potential are considered, sub-Planckian values of the inflaton field are required and 
adequate reheating of the universe is accomplished via curvature-induced ll24l couplings of the inflaton 
to matter fields. Comparisons with the results obtained for the minimal version of both inflationary 
models are also displayed. 

Below, we describe the generic formulation of non-MI (Sec. [2]) and then apply the relevant results, 
for appropriate choices of f(a), in the case of non-MCI and non-MHI in Sec. [3] and [4] respectively. 
Finally, Sec. [5] summarizes our conclusions. Throughout the text, we set natural units for the Planck's 
constant, Boltzmann's constant and the velocity of light (h = c = &b = 1) the subscript , x denotes 
derivation with respect to (w.r.t.) the field \ ( e -g-> ,xx = d 2 /dx 2 ) and a bar over a field x denotes 
normalization w.r.t the reduced Planck mass, mp = 2.44 • 10 18 GeV, i.e., x = x/ m P- Finally, we 
follow the conventions of Ref. |[25l for the quantities related to the gravitational sector of our set-up. 

2 Inflation with non-Minimal Gravitational Coupling 

Non-MI, by its definition, can be realized by a scalar field non-minimally coupled to Ricci scalar 
curvature. The formulation of a such theory is described in Sec. 12.11 Based on it, we then derive the 
inflationary observables and impose observational constraints in Sec. |2.2j 
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2.1 NON-MINIMALLY CURVATURE-COUPLED SCALAR THEORY 

The dynamics of a scalar field a non-minimally coupled to 1Z through a coupling function f(a) is 
controlled, in the Jordan frame, by the following action - see, e.g., Ref. |4): 

S = j d^x^g L± m \f{a)n + \g^d^d v a - V {a)\ , (2.1) 

where g is the determinant of the background Friedmann-Robertson- Walker metric [25]. To guarantee 
the validity of the ordinary Einstein gravity at low energy, we require /((c)) = 1, where (a) is the 
vacuum expectation value (v.e.v) of a at the end of non-MI. 

The action in Eq. (12.11 ) can be brought in a simpler form by performing a conformal transformation 
ll26l to the so-called Einstein frame where the the gravitational sector of our model becomes minimal. 
Indeed, if we define the Einstein-frame metric 



-g = f 2 V=g and r u = gT/f, 

n={n + 2,U\uf + 3g^d„fd„f/2n /f 



- _ f j y y — j v y a " u y — y i j i a n\ 

g^v — j g^v \£ _ lrn , orn ,„ j. , o„uvx .fa t ir>{2\ it y --- ] 



- where □ = \—g) (yf—gd^) and hat is used to denote quantities defined in the Einstein frame - 

and introduce the Einstein-frame canonically normalized field, a, and potential, V, defined as follows: 

the action in Eq. (12.11 ) can be simplified, taking the form 

S = f d 4 x^ (-~mlK + \r v d^d u a - V (a)\ . (2.4) 

Based on the action above, we can proceed readily to the analysis of non-MI in the Einstein frame 
using the standard slow-roll approximation (9j|27l - see below. It can be shown ll28l that the results 
calculated this way are the same as if we had calculated them with the non-minimally coupled scalar 
field in the Jordan frame. 

One of the outstanding features of the scalar theories with non-minimal f(a) is that a can decay via 
gravitational effects ll24l even without explicit couplings between a and matter fields. This is, because 
couplings arise spontaneously when a settles in its v.e.v, (a), and oscillates, with coupling constants 
involving derivatives of f(a) calculated for a = (a). If we identify a as the inflaton, these couplings 
can ensure the reheating of the universe. Assuming the existence of a bosonic field minimally coupled 
to gravity, with negligible mass compared to the mass of a, m a = V )aa ((a)) l l 2 , we get M241I29I for the 
reheat temperature 

is the decay rate of a, in the regime T <C m a which is valid in our applications. Clearly, this construc- 
tion is applicable if / j(T ((a)) ^ (and this is valid for the /(cr)'s considered in Sec. [3] and [4]). Also, 
assuming the particle spectrum of SM, we set g p * = 106.75 for the relativistic degrees of freedom. 

2.2 Inflationary Observables - Constraints 

Under the assumption that (i) the curvature perturbations generated by a is solely responsible for the 
observed curvature perturbation and (ii) there is a conventional cosmological evolution (see below) 
after inflation, the inflationary parameters can be restricted imposing the following requirements: 
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(a) The power spectrum P-ji of the curvature perturbations generated by a at the pivot scale 
K = 0.002/Mpc is to be confronted with the WMAP5 data H21, 

pV°- 1 ^) 3/2 - l J ^)l ^^ lO! 1Q-5 (26) 
n 2V3™3 \Vf(d*)\ 2V3vrm3 \V a (a,)\ 

where er* [a*] is the value of a [a] when k* crosses outside the inflationary horizon. 

(b) The number of e-foldings, N* , that the scale fe* suffers during FHI is to account for the 
total number of e-foldings N tot required for solving the horizon and flatness problems of standard big 
bag cosmology, i.e., TV* = iV to t. Specifically, we calculate iV* through the relation 

i r* ^ v i r* 9 V 

N*= ^ da^ = ^ da J 2 ^, (2.7) 

where <7f [d{] is the value of a [a] at the end of inflation, which can be found, in the slow-roll approxi- 
mation and for the considered in this paper models, from the condition 

max{?(cjf), |?7(cf)|} = 1, where 

ml ( V\aY nil /y„\ 2 2 Vffff ml ( V aa V ff J CT \ 

and r] = m P = — f -k — - ■ (2.8) 

v J 2 \ V v J ) 

The required A^ tot at can be easily derived |[T9l consistently with our assumption of a conventional 
post-inflationary evolution. In particular, we assume that inflation is followed successively by the 
following three epochs: (i) the decaying-inflaton dominated era which lasts at a reheat temperature T r h, 
(ii) a radiation dominated epoch, with initial temperature T r h, which terminates at the matter-radiation 
equality, (iii) the matter dominated era until today. In particular, we obtain - c.f. Ref. |fl9ll 

where the last term emerges from the transition from the Jordan to Einstein frame. Note that 
R = yffR with R being the scale factor of the universe. 

(c) The (scalar) spectral index, n s , its running, a s , and the scalar-to-tensor ratio r are to be 
consistent with the fitting (HI of the WMAP5 plus BAO and SN data, i.e., 

(a) n s = 0.96 ± 0.026, (b) - 0.068 < a s < 0.012 and (c) r < 0.22, (2.10) 

at 95% confidence level (c.L). The observable quantities above can be estimated through the relations: 

2 ^ 
n s = 1-6?* + 277*, a s = - (4^ _ (n s - l) 2 ) - 2$, and r = 16e, (2.11) 

where £ = m\V jV 2 = rripV^rj^/V J 2 + 2ffe and the variables with subscript * are eval- 
uated at a = a*. Note, in passing, that the utilized here non minimal /(cr)'s do not produce |[30l 
observationally interesting non-gaussianity - for reviews see, e.g., Ref. [31 ]. 

(d) To avoid corrections from quantum gravity, we impose two additional theoretical con- 
straints on our models - keeping in mind that V(af) < V{a*)\ 

(a) F(cr*) 1/4 < m P and (b) a* < m P . (2.12) 

Although it is argued |[20l[32l that violation of Eq. (12.12b ) may not be necessarily fatal, we insist on 
imposing this condition in order to deliberate our proposal from our ignorance about the Planck-scale 
physics. To be even more conservative, we have to check the hierarchy between the ultraviolet cut-off, 
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A, of the effective theory and the inflationary scale. The former can be found from the non-renormali- 
zable terms arising in Eq. ( 12.41 ), whereas the latter is represented by V{a*) l l 4 or, less restrictively, 
by the corresponding Hubble parameter, = V "(a*) 1 / 2 / 'y/3mp. In particular, the validity of the 
effective theory implies @ 

(a) y(<7*) 1/4 < A or (b) H* < A. (2.13) 

This requirement applies mainly in cases where the involved in /(<?) constant c-jz takes relatively large 
values - as for SM non-MI [2]44j - jeopardizing, thereby, the validity of the classical approximation, 
on which the analysis of the inflationary behavior in this section is based. 

3 Non-Minimal Chaotic Inflation 

We focus on CI driven primarily by a quadratic potential of the form 

V = \m 2 a 2 + V [C where V rc = In ^ (3.1) 

are. radiative corrections l23l to the inflationary potential. The bulk of our results - see Sec. 13.21 - are 
independent of the renormalization scale, Q, which is set equal to rop. We below recall (Sec. 13.11 ) 
the results for MCI (with f(a) = 1) and describe (Sec. 13.21 ) our findings for non-MCI, adopting the 
following coupling function - recall that a = a/mp: 

f{p) = (1 + c n a)~ n with n = ±1. (3.2) 

Note, in passing, that results for non-MCI with quartic potential (V = Ao" 4 /4!) are presented 
in Ref. fTl 1551 1531 . The inflationary scenario based on this potential with f(a) = 1 seems to be 
excluded [12,20] due to the enhanced predicted r. As we explicitly verified, if we employ the standard 
non-minimal coupling function, f(a) = 1 + c-jia 2 , with 80 < c-ji < 300 and 0.2 < A/10 -4 < 3 - c.f. 
Ref. B331I341 - we can rescue the model consistently with the constraints of Sec. l2.2l for an indicative 
T r h = 10 10 GeV. In particular the lower [upper] bound of the allowed regions of c-ji and A comes from 
Eq. ( 12. 12b ) [Eq. (12. 13a ) with A = mp/cyj]. Note, however, that the standard non-trivial f(a) does not 
support reheating along the lines of Eq. (I2.5I ). 

3.1 Results for MCI 

For MCI the slow-roll parameters and the number of e-foldings suffered from k* can be calculated 
applying Eq. (12.81 ) and Eq. (12.71 ) - after removing hats and setting J = 1 - with results 

e = r] = 2/a 2 and N* = (a 2 - a 2 ) /4. (3.3) 

Using these results, imposing the condition of Eq. (12.81 ) and employing Eq. (12.1 11 ) we can derive 

CT f = V2, a* ~ 2 v / iV*, n s ~ 1 - 2/iV* and r ~ (3.4) 

Clearly trans-Planckian values of a are required and observationally favored n s and r are obtained. 
More precisely, imposing the requirements (a) and (b) of Sec. l2.2l for several T r h's we get numerically 
the values of <r*, m, n s , a s and r listed in Table [T]- c.f. Ref. |[20l . As T r h decreases, iV* decreases too 
- see Eq. ( 12.91 ) - and so, a* and n s slightly decrease whereas r increases - see Eq. (I3.4I ). The resulting 
n s , a s and r lie within the range of Eq. ( I2.10I ). In all cases, Eq. ( 12. 12a ) is valid whereas the upper 
bound of Eq. ( 12. 12b ) is surpassed. 
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T rh (GeV) 


cr*/mp 


m (10 13 GeV) 


n s 


a s (10- 4 ) 


r 


10 10 


15.13 


1.6 


0.965 


6.1 


0.139 


10 6 


14.73 


1.69 


0.963 


6.7 


0.147 


10 5 


14.61 


1.72 


0.962 


7 


0.15 


10 4 


14.5 


1.74 


0.962 


7.2 


0.152 



Table 1: Values of parameters allowed by Eqs. H2.6T) . (2.7\) and (2.9b for MCI with several T r h 's. 



3.2 Results for non-MCI 



From Eqs. (12.81) . (13.31 ) and (13.41 ) we can infer that the amplitude of the infiaton field within non-MCI 
can become sub-Planckian if J ~ l//(er) S> 1 and V a /V ~ V a /V. These two objectives can be 
achieved if we employ f(o~) given by Eq. ( 13.21 ) with n > and c-r 3> 1. Another possibility would 
be to take f(a) = exp (— oro) with or ~ 10. However, in the latter case the resulting r violates 
Eq. ( 12. 10b ) and therefore, this option can be declined. Similar problem arises also if we use n > 1 
- see Sec. 13.2.11 On the other hand, for n = —1, V in Eq. (12.31 ) becomes very flat for sufficiently 
large er's and so, a new type of non-MCI can takes place. Decreasing n for n < we find inflationary 
solutions, only for or < 0.001, which break Eq. (12.12b ). Similar conclusions are also drawn for the 
standard non-minimal f(a) - see Ref. 11331 . 

In our numerical code we use as input parameters m, a*, or and n. For every chosen n and c-r 
we restrict m and cr* so as the conditions (a) and (b) of Sec. 12.21 - with T r h evaluated consistently with 
Eq. ( 12.51 ) - are fulfilled. Our results for n = +1 and n = —1 are presented respectively in Fig. Q] and 
Fig. [2l where we draw the allowed values of m (solid line) and T r h (dashed line) [erf (solid line) and cr* 
(dashed line)] versus c-r for non-MCI (a) [(b)]. For both n = ±1, satisfying Eq. (12.12b ) gives a lower 
bound on or - see Fig. [B(b) and Fig. |2]-(b). On the other hand, the upper bound on c-r comes from 
Eq. (12.10b ) for n = +1 and from the fact that the enhanced resulting m's destabilize the inflationary 
path through the radiative corrections in Eq. (13.11 ) for n = — 1. From our data we also remark that the 
resulting m's are almost two orders of magnitude lower [larger] than those obtained within MCI for 
n = +1 [n = — 1]. These results depend, though very weakly, on A^ot and therefore, on the reheating 
mechanism - see Eq. (12.9I ). All in all, we obtain 

625 <c n < 2.1 • 10 7 , 47 > — > 1.6 and 52 > N* > 47.9, for n = -1, (3.5) 

771 

83 < c n < 3120, 3 < 1Ql2 ^ < 8.6 and 58.8 < N* < 59.9, for n = +1. (3.6) 

In both cases, the predicted n s and r lie within the allowed ranges of Eq. (12.10b ) and Eq. (12.10b ) 
respectively, whereas a s remains quite small. Our numerical results can be interpreted through some 
simple analytical expressions which are presented in Sec. 13.2. T\ TSec. [3T2.2H for n = +1 [n = — 1]. 
There, we also comment on the naturalness of our models, following the arguments of Ref. (6j|71. 

3.2.1 Non-MCI with n = +1 

To justify our choice for the negative exponent in Eq. ( 13.21 ) we present our formulae below for a general 
n > 0. Substituting Eq. ( 13.21 ) into Eqs. ( 12.31 ) and ( 12.51 ) and taking into account that c-r » 1, we obtain 

1 ooA a \ 2n m 2 cl n a 2 ^ .„ 1 ml 

(3.7) 



J~^, F = - m V 1 + CK- ~ » and T 

v 2 \ mp / 2mp 



192?r ml 
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(a) (b) 

Figure 1: The allowed by Eqs. (2.5\) , (12.6D . (2.7\) and (2.9\) values of m (solid line) and T r h (dashed 
line) [erf (solid line) and er» (dashed line)] versus c-r for non-MCI with n = — 1 (a) [(b)]. The gray 
segments denote values of the various quantities fulfilling Eq. i2.13b ) too. 



where m a = m and obviously (<r) = 0. Upon use of Eqs. (12.81 ). (12.71 ) and (13.71 ). the slow roll parameters 
and iV* read 

^ 2(1 + n) 2 ^ 2(l + n)(l + 2n) (l + 2nU c ^ (a^+ 2 - a^ +2 ) 

(a) ^^ \ ' , ^ ^ = \^ > e and (b) N* ~ | * f / • (3.8) 

C ^, CJ d^a n+j! (1+n) 2(l+n)(2 + n) 



Imposing the condition of Eq. (12.81 ) and solving then Eq. (13.8b ) w.r.t a* we arrive at 

a f ~ (2(l + 2n)(l + n)/c^) 1/( " +2) and ct* ~ (2(1 + n)(2 + n)iV,/c^) 1/( " +2) • (3.9) 



Inserting the last results into Eq. (13.8b ). we find through Eq. (12.111 ) 

(a) n s ~ l-3e* = 1-3(1 + n)/(2 + n)iV* and (b) r ~ 16(1 + n)/(2 + ra)JV*. (3.10) 

Letting vary within its allowed region for n = +1 - see Fig.Q]- we find ra s ~ (0.952 — 0.955) and 
r ~ (0.2 — 0.22). Clearly, increasing n leads r above the range of Eq. (12.10b ). Therefore, we hereafter 
concentrate on n = +1 which assures an observationally safe and, at the same time, exciting r. 

Comparing our findings with those obtained for MCI - see Table Q]- we notice that the resulting 
here n s 's are a little lower, whereas r is significantly elevated and can be probed in the near future from 
the measurements of PLANCK satellite |[35ll . Note, in passing, that the so-called Lyth bound |f36l on 
the ct variation, Act, gets modified within non-MI. Namely, combining Eqs. ( 12.71 ) and ( 12.81 ) we find 

da prmp . PFmp A ~ . , — top 

— = J-—f => Act = A /--f AiV Act > V27^- ~ (2.5 - 0.083) top/100, (3.11) 
eJjV y 8 J \ 8 J J 

taking ||36l AiV ~ A A = 4 and assuming negligible variation of /(ct) from its value at ct = ct*. 
Therefore, large r's do not correlate necessarily with trans-Planckian Act's within non-MI. On the 
other hand, ct as evaluated from Eq. ( 12.31 ), ct ~ \J ciza 3 /mp, remains trans-Planckian. 

The resulting V in Eq. (13.71 ) is non-renormalizable and suggests that the theory breaks down for 
energies of the order A = mp/c^. Checking the consistency with Eq. (12.13b ) we find numerically: 

0.03 < H*/A < 1 for 625 < c n < 2.26 • 10 4 and 1 > ct* > 0.3, (3.12) 
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10 2 r 



E 10: 




1 4±J 



K 

(a) 



K 

(b) 



Figure 2: The same as in Fig.{l\but forn = — 1. The gray segments here denote values of the various 
quantities fulhlling Eq. (2.13k ) too. 



where the corresponding ranges of values are depicted by the gray segments of the lines in Fig. [TJ The 
range in Eq. (13.121 ) turns out to be a little more comfortable than the one we get within SM non-MI - 
c.f. Ref. @. However, Eq. (IXT3"k ) is violated, since y(cr*) 1 / 4 /A > 5.8. 

On the other hand, non-renormalizable terms in the action of Eq. ( 12. 11 ) and (12.41 ) indicate that 
A = mp. In fact, such terms arise from the first term in Eq. (12.11) and the second one in Eq. (12.3b - The 
form of these terms is generated expanding the relevant coefficients in series around o = with the 
following result - an expansion in the small field limit, c-jia <C 1, fails to reproduce the exact results: 

1-3^ + 10 f^Y + - N 



m%flZ 



mp 



c-r '* 



(3.13a) 



fl 

and mp-^g^df+o-dvo- B 



1 a la 

-tt 1-8— + 45 — 



+ ---\srdaad v a, (3.13b) 



where h^ v denotes the graviton field involved in the expansion [6,7] of the metric g^ v ~ il^ + h^/mp 
around the Minkowski space with metric t]^ u and 1Z is approximated linearly. Given these ambiguities, 
we do not consider Eq. (I3.12I ) as absolute constraint. 

3.2.2 NON-MCI WITH n = — 1 

A completely different situation from that studied in Sec. 13.2.11 emerges for n = —1 in Eq. ( I3.2I ). 
Indeed, substituting Eq. (13.21 ) into Eqs. ( 12.31 ) and (12.51 ) and taking into account that c-ji 3> 1, we obtain 

J ~ ^/372 a-\ V ~ and T a ~ (3-14) 



2 2 

m nip 

— ^ and T 



192vr mp 



where m a = m and obviously (a) = 0. We observe that V exhibits a flat plateau as we obtain for the 
quatric potential with the standard non-minimal f(a) - c.f. Ref. |2]-I4l[34l. Employing Eqs. (12. 8b . (12.71 ) 
and (13. 14b . the slow roll parameters and A^* read 

4 4 . ^ .... ft 3cr, 



a e ~ 



-tcfia » — e and (b) 



(3.15) 



As opposed to our findings in Eqs. (13.31 ) and (13. 8b . notice that r) < here. Imposing the condition of 
Eq. (12.8b and solving then Eq. (13.15b ) w.r.t we arrive at 



<7f ~ 2/V3oji and a* ~ 4 N*/3cn- 



(3.16) 
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Inserting the last results into Eq. ( 13.15k ), we find through Eq. ( 12.111 ) 

(a) n s ~ 1 + 2% = 1 - 2/iV* ~ (0.967 - 0.97) and (b) r ~ 12/A\ 2 ~ (0.002 - 0.003), (3.17) 

where the ranges above are derived numerically letting c-ji vary within its allowed region - see Fig. |2j 
Notice that the resulting n s 's and r's are identical to those derived in Ref. EH. Comparing them with 
those listed in Table Q] or given in the paragraph below Eq. (13.101 ) we remark that r is significantly 
reduced, whereas n s is close to the value obtained in MCI and a bit larger than the one extracted for 
non-MCI with n = +1. 

As for the latter case, non-renormalizable terms in the action of Eq. (12.11 ) indicate an effective 
cutoff A = mp/c-jz, since 

mlfn [afh^erhl + a^a^ + a^h) o.is) 

with h = h 1 ^ = h^. On the other hand, the second term in Eq. (12.31) gives exactly the same result as 
in Eq. (13. 13bb since f :(7 /f is identical for both n = ±1 in Eq. (13.2I ). Checking the consistency with 
Eq. ( 12. 13a ) we find numerically 

0-3 < V(a.)V*/A < 1 for 83 < c n < 313 and 1 > <r» > 0.27, (3.19) 

where the corresponding ranges of values are depicted by the gray segments of the lines in Fig. [2] On 
the other hand, Eq. ( |2.13b ) is satisfied in the whole parameter space of these figures. Consequently, 
non-MCI with n = — 1 can be characterized as more natural than the one with n = +1. 

Concluding this section, let us emphasize that, in contrast to the models suggested in 
Ref. [20], non-MCI is not of hilltop type and so, complications related to the initial conditions are 
avoided. Indeed, for non-MCI with n = +1,V remains concave upwards, whereas with n = — 1, V 
develops a plateau without distinguished maximum. As we explicitly checked, possible inclusion of 
extra radiative corrections in Eq. (13.11 ) due to a coupling of a to fermions - considered in Ref. ll20l - 
do not affect our proposal for values of the relevant Yukawa coupling constant, h, lower than about 
10 -3 . For such /i's, the decay width of the inflaton due to this channel dominates over the one given 
by Eq. (1231 



4 Non-Minimal Hybrid Inflation 

Hybrid inflation can be realized in the presence of two real scalar fields, a and <f>, involved in the 
following potential ifTTl 

V(^a) = K 2 (M 2 -^) +—- + ^L—, (4.1) 

where M, m are mass parameters and k, A are dimensionless coupling constants. The global minima of 
V lie at ((a), (</>)) = (0, ±2M). Therefore, V leads to a spontaneous symmetry breaking of a global 
or local symmetry depending on the nature of the waterfall field <j). In the latter case, topological 
defects may be also produced via the Kibble mechanism [37]. Trying to keep our approach as simple 
as possible we below assume that this is not the case. 

In addition, V in Eq. (14.11 ) gives rise to HI. This is because V possesses an almost c-flat direction at 
= with constant potential density equal to Vq = V (4> = 0, a) = k 2 M 4 , for m = 0. The effective 
mass squared of the field <fi along this direction is 



ml = -k 2 M 2 + AV/2 > <S> a > a c = \/2kM/\ . 



(4.2) 
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Thus, for o > o c the = direction represents a valley of minima which can serve as inflationary 
trajectory. On this path the potential of HI takes the form 

V = V + X -m 2 o 2 + V IC (4.3) 
where V rc is the one-loop correction (to the tree-level potential) which can be written as 

1 ( .4, ™ 2 , J2rr f„ ^ K 2 M 2 , ^\ ... ( 1 



V IC = ^ [rn 4 ln-^ + K 2 V (x-irin-^(x-l)j with x = l—j (4.4a) 

n 2 V ( k 2 M 2 3\ 

' x 2 In — =- + - for x > 1. (4.4b) 



64vr 2 V Q 2 2 

Here, Q is a renormalization scale which can be conveniently chosen ||2T1 equal to a c which practically 
coincides with the value of a at the end of HI, at, for both MHI and non-MHI. 

We below review (Sec. 14.11 ) the results for MHI (with f(a) = 1) and describe (Sec. 14.21 ) our 
findings for non-MHI, seeking the following non-minimal coupling function for the inflaton - for 
earlier attempts on non-MHI, see Ref . 11381 : 

f(a) = l-c n a/(l + a) 2 (4.5) 

where we use, as usually, the shorthand a = a/mp. As regards the waterfall field we can assume that it 
is either minimally coupled to gravity or its coupling function is f(<p) since /(0) = 1 and f(2M) ~ 1 
for cji <C 1 and M < mp. 

4.1 Results for MHI 

We can get an impression of the expected results for MHI, if we calculate the, involved in the inflation- 
ary dynamics, derivatives of V in Eq. (14.3I ). Namely we have 



V a = m 2 a + ^(x-l)(2ln^(x-l) + l). (4.6) 



32vr 2 cr v ; V Q 

We observe that there are two contributions in V a . The first one arises from the tree-level potential 
whereas the second one comes from the radiative corrections in Eq. (I4.4ab . When the first contribution 
dominates over the second one, we obtain the well-known tree-level ifTTTl results, 2Vt r * and rj tr , for iV* 
and rj respectively - note that we identify <jf with a c : 

1 Gif 2 m2 

N tr * = — m — with rjtr = mp— > e . (4.7) 

In this regime, the resulting n s clearly - see Eq. (12.1 II ) - exceeds slightly unity in contrast to the 
observationally favored results of Eq. (12. 10b ). Moreover, as we find numerically, the lower k and/or m 
we use, the closer <r* is set to a c - This is the first kind of tuning occurred within MHI. 

Nonetheless, taking into account that the logarithm in Eq. (14.61 ) turns out to be negative, we can 
show that, for every m, there is k such that V develops a maximum at a = <r max , which can be 
estimated by numerically solving the condition V a (a max ) = 0. At a = cr max , V a(T given by 

becomes negative and so, r] and n s start decreasing for <t* close <J max - see Eqs. ( 12.81 ) and (12 -lib . As 
for any model of hilltop inflation, the lower n s we obtain, the closer cr* is located to <7 max . This is 
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Figure 3: The allowed by Eqs. $2J$, HZ7\) and %2J% values of n s (a) and M (b) versus k for MHI 
with T r h = 10 10 GeV, n = A and several m's indicated in the graphs. Shown are also in the table the 
allowed ranges of the various parameters for n s in the range ofEq. i2.10k ) limited by thin lines (a). 



a second kind of tuning which remains even for non-MHI - see Sec. 14.21 To quantify somehow the 
amount of the tunings encountered in the considered model, we define the quantities: 

(a) A m * = and (b) A c * = (4.9) 

The above rough estimations can be verified by our numerical computations. In our code, we use 
as input parameters k, A, m, M, a* and T r h- In our analysis for MHI, we fix T r h = 10 10 GeV and 
k = X - possible variation of these two choices do not modify our conclusions in any essential way. 
For any chosen k and m we then restrict M and cr* so as the restrictions (a), (b) and (d) of Sec. 12.21 
and Eq. ( 14.21 ) are fulfilled. Using Eq. (12.1 11 ) we can extract n s , a s and r. Our results are presented in 
Fig- E]-(a) [Fig. [3]-(b)] where we design the allowed values of n s [M] versus k for m = 1 TeV (solid 
line) or m = 10 3 TeV (dashed line) or m = 10 6 TeV (dot-dashed line) or m = 10 9 TeV (dotted line). 
The region of Eq. (12.10a ) is also limited by thin lines. The various lines terminate at low k's due to the 
saturation of Eq. (12.12b ) and at large k's since the imposed conditions can not be fulfilled. 

Clearly, the almost horizontal part of the various lines, which exceeds the observational limits of 
Eq. (12. 10a ), in the k — n s plane corresponds to the dominance of the tree-level potential. However, for 
any m and relatively large k's we can obtain acceptable n s 's even without inclusion of extra radiative 
corrections due to a possible coupling of the inflaton to fermions - c.f. Ref. (2TJ. On the other hand, 
it is worth emphasizing that the allowed range of k's for each m is severely tuned. Indeed, confining 
n s within the range of Eq. (12.10a ) we find the ranges of the parameters listed in the table of Fig. [3] 
From the outputs there, we also remark that k's, M's, A^'s and A m *'s increase with m. Therefore, 
the natural realization of MHI requires large m's. In this case too, M turns out to be well above its 
value within the SUSY version HI - c.f. Ref. II15N16I . Needless to say, finally, that the resulting a s 's 
and r's turn out to be vanishingly small and so, uninteresting. In conclusion, MHI (with the minimal 
possible radiative corrections) is rather disfavored by the current observational data. 
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4.2 Results for non-MHI 

From the analysis of MHI we can deduce that reduction of n s for a wider range of k's can be achieved 
if the slope of V becomes steeper. This objective can be achieved if we employ f(u) given by Eq. (14.51) 
with en <C 1. Another possibility would be /(<j) = exp (— c^cr) or that of Eq. (13.21) with n > and 
c-ji ~ 0.1. However, in these cases the resulting cr* violates the bound of Eq. (12. 12b ) and therefore, 
these options are not adoptable. Moreover, imposing on non-MHI with the standard non-minimal f(a) 
the constraints (a), (b) and (d) of Sec. |2.2| and Eq. ( 14.21 ), we are obliged to use a tiny c-ji ~ — 10 -3 , which 
has no sizable impact on reducing n s . Consequently, this last choice can not become observationally 
viable, too. 

Differentiating Eq. (14.51 ) w.r.t a, substituting into Eqs. (12.31 ) and (12.51 ) and taking into account that 
c-ji <C 1, we obtain 

/„ = + % f„ = 2 J {2 ~% V ~V , J ~ 1 and T CT = ^4 (4-10) 

' m P (l + a) 3 ' to 2 (1 + ct) 4 128vr to 2 



where m a = \/2X 2 M 2 + to 2 . Despite the fact that V given by Eqs. (12.31) and (14.31) is practically equal 
to Vq - since f(a) ~ 1 for a *C 1 -, its inclination is mostly dominated by the term — 2Vo/ j(T of V >a . 
Indeed, upon use of Eqs. (f278b . (I4.4bb and (14.101 ) we find 

£ = T 1 (- 2/ - + -^ + l6A 1 "^j • <411) 

In a sizable portion of the parameter space, the first contribution to e in Eq. (14.1 II) overshadows the 
others two. As a consequence, V develops a maximum at ct = cr max for (<7 max ) = 44> a max ~ 1 
with V aa (cr max ) < 0. In fact, inserting Eqs. (12.31 ) and (14.31 ) into Eq. (12.81 ) we end up with 

n / m 2 3k 2 x 2 , k 2 M 2 \ 

, = - 2 ^-2/^ + ^ + ^ In J , (4.12) 

which is negative for dominant / jCTCT with ct < 2. Combining Eqs. (14.121 ) with (12.11a ) we can easily 
infer that c-ji > for a < cr max strengthens significantly the reduction of n B . Neglecting the two last 
terms in the right-hand side of Eq. (14. lib , we can estimate iV* via Eq. (12.71 ) with result 

N* — 7T—2 [ C ~7~ = ~F~ ( ( 21 + + ^* " ( 21 + 6 ^c + ^c) + 241n ^— ^ ] • (4.13) 

As we verify numerically, the formula above gives accurate results for m < 10 6 TeV and sufficiently 
low k's. However, since a* depends on N* in a rather complicate way, it is not doable to find an 
analytical result for n s as a function of iV* - c.f. Eq. (13.91 ) and Eq. ( I3.9I ). Therefore, our last resort is 
the numerical computation, whose the results are presented in the following. 

In our code, we use as input parameters k, A, to, M, a* and c^. Note that T r h is calculated via 
Eq. ( |2.5b . For every chosen k, A, to and c-ji, we can restrict M and cr* so as the conditions (a), (b) 
and (d) of Sec. 12.21 and Eq. ( 14.21 ) are fulfilled. Through Eq. ( 12.111 ) we can then extract n s and a s . 
Following this strategy, in Fig. |4]-(a) [Fig. |4]-(b)] we display the allowed values of n s [M] versus or. 
with to < 10 6 TeV, k = 10 -5 and A = k (solid lines) A = 5k (dashed lines) and A = 0.5k (dotted 
lines). The region of Eq. (12.10a ) is also limited by thin lines. We observe that as c-ji increases, n s 
decreases entering the observationally favored region of Eq. (12.10a ). On the other hand, M increases 
with c-jz until a certain oji ~ 0.03 — 0.05 and then decreases. Surprisingly the value of on, at which 
the maximum M is encountered, corresponds more or less to the central observational n s ~ 0.96. 
We also observe that increasing A above k with fixed on, n s drops but M raises. These results can be 



13 



Non-Minimally Gravity-Coupled Inflationary Models 




Figure 4: The allowed by Eqs. JZ6b . fZ7\) and - with T xh given by Eq. iO - values of n s 
(a) and M (b) versus c-r. for non-MHI with m < 10 6 TeV, k = 10~ 5 and several A's indicated in the 
graphs. The region ofEq. (2.10k ) is also limited by thin lines. 



understood as follows: As A/re elevates a c decreases - see Eq. (14.21 ) - and therefore, cr* decreases, with 
fixed iV*. This effect causes an increase of \f,a(cr*)\ and f )ffa ((T*) - see Eq. (14.101) . As a consequence, 
M increases too, since M is proportional to fj due to Eq. (12.6b - Also, |r/| increases - according to 
Eq. (14.121) - and so, n s drops efficiently - see Eq. (12- lib . 

Confronting non-MHI with all the constraints of Sec. 12.21 consistently with Eq. (14.21 ). we can de- 
lineate the allowed (lightly gray shaded) regions in the k — cr [re — M] plane as in Fig. (5]-(ai), (bi) 
and (ci) [Fig. [5j-(a 2 ), (b 2 ) and (c 2 )]. In Fig. |5]-(ai) and (a 2 ) we take A = re. Our results for this choice 
are m-independent for any re and m < 10 6 TeV. On the other hand, in Fig.[5]-(bi) and (b 2 ) [Fig. f5]-(ci) 
and (c 2 )] we set A = 5k and m = 10 s TeV and [A = 0.5k and m = 10 7 TeV]. The conventions 
adopted for the various lines are also shown in the left-hand side of each graph. In particular, the gray 
dot-dashed [dashed] lines correspond to n s = 0.986 [n s = 0.934], whereas the gray solid lines have 
been obtained by fixing n s = 0.96 - see Eq. (12. 1 Ob - For re's below the solid black line, our initial 
requirement in Eq. (12.12b ) is violated. For re's larger than those depicted in the graphs we do not find 
solutions consistent with the imposed restrictions of Sec. 12.21 The upper bounds of the allowed regions 
in the k — M plane come from or leading to n s = 0.96 - see Fig- SJ-(b). Although this result may not 
rigorously correct, it is accurate enough for our pictorial purposes. In all cases, the allowed ranges of 
k's - although restricted to values lower than 0.001 - are much more wide and natural than the ones 
obtained for MHI - c.f. Table of Fig.[3j Confining n s to its central observational value, we obtain the 
ranges of the various parameters arranged in the Table of Fig. [5] We observe there that, for fixed ra s and 
increasing k, c-ji and M decrease whereas A c * and A m * increase. As a consequence, for any m, the 
tuning regarding A C! „ is greatly alleviated compared to the outputs of MHI, whereas we are let with the 
usual mild tuning required for A m *. This is present to any inflationary hilltop model - c.f. Ref. [ 16]. 
The allowed M's mostly exceed the SUSY grand unification scale, Mqut — 2.86 • 10 16 GeV, whereas 
T r h mostly increases with k, as can be noticed via Eqs. ( 12.51 ) and d4- 10b - 

From our findings, we can conclude that: (i) the required c^'s are rather low and so, complications 
related to the hierarchy between the inflationary scale and the effective cutoff of the theory are avoided; 
(ii) our results depend rather weakly on the variation of m, for m < 5 • 10 8 TeV; (iii) as m raises above 
5 • 10 8 TeV and k drops below 0.001, Eq. ( 12. 12b ) is eventually violated and so, our scheme becomes 
unapplicable; (iv) similarly to MHI, a s and r turn out to be negligibly small. 
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Figure 5: Allowed (lightly gray shaded) by the restrictions of Sec. \2.2\ consistently with Eq. (4.2\) 
areas in the k — or. [k — M] plane (an,, bi and ci) f(a 2 b 2 and c 2 )J for non-MHI. We take n = A and 
m < 10 6 TeV (ai and a 2 ) orm = 10 8 TeV and A = 5k (bi and b 2 ) orm — 10 7 TeV and A = 0.5k 
(bi and c 2 ). The conventions adopted for the various lines are also shown. The allowed ranges of the 
various parameters forn s = 0.96 are listed in the table. 
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As in the case of non-MCI, our proposal remains intact even if we add fermion-dominated one- 
loop radiative collections in Eq. ( 14.31 ) - c.f. Ref. EH - provided the values of the relevant Yukawa 
coupling constant, h, remains lower than about 10~ 4 . For h's close to this value, the decay width of 
the inflaton, due to this channel dominates over the one given by Eq. (12.5I ). 

5 Conclusions 

We considered the non-SUSY version of CI (driven by quadratic potential) and HI, assuming a non- 
minimal coupling function, f(a), between the inflaton field and the Ricci scalar curvature. Using 
the freedom of choosing this scalar function, we deliberated CI from the problem of trans-Planckian 
inflaton values and achieved observationally acceptable n s 's for a wide range of the parameters of HI. 
As a bonus, the selected /(c)'s give rise to Yukawa-type interactions between the inflaton and matter 
fields leading to a successful post-inflationary reheating. Afterwards, the proposed /(<r)'s reduce to 
unity and so, the Einstein gravity is naturally recovered. 

Specifically, the adopted forms of f(a) are given by Eq. (13.21 ) and Eq. (14.51 ) for non-MCI and 
non-MHI, respectively. In both cases, the parameter c-jz involved in f(a) can be constrained so as 
the results of the inflationary models can be reconciled with a number of theoretical and observational 
restrictions. Our results are as follows: 

• In the case of non-MCI, we find 625 < c n < 2. MO 7 resulting to n s ~ 0.955 and r ~ (0.2-0.22) 
for n = +1 and 83 < c n < 3120 resulting to n s ~ 0.967 and r ~ (0.002 - 0.003) for n = -1. 
In sharp contrast to MCI, only sub-Planckian values of the inflaton field in the Jordan frame are 
utilized avoiding, thereby, destabilization of the inflationary scenarios from possible corrections 
caused by quantum gravity. Comments on the naturalness of the models are also given. 

• In the case of non-MHI, the chosen f(a) leads to hilltop-type inflation for a wide range for k's. 
As a consequence, observationally acceptable results require a proximity between the values of 
the inflaton field at the maximum of the potential and at the horizon crossing of the pivot scale. 
The amount of this tuning was measured by the quantity A mt defined in Eq. ( 14.9b ). E.g., for 
m < 10 6 TeV and the observationally central value of n s , we find cn c± (0.015 — 0.078) with 
M ~ (1 - 16.8) • 10 17 GeV, A = k ~ (2 • 10~ 7 - 0.001) and A m * ~ (0.91 - 32)%. Compared to 
MHI, we find that the observational requirements can be satisfied without tuning severely neither 
k nor A c * defined in Eq. ( 14.9k ) even for low m's - see Tables of Figs. [3] and [5j Therefore, the 
proposed non-MHI is more favored by the current data. 

We explicitly checked that, for both models of non-MI, the proposed scheme remains valid even 
if an extra coupling of the inflaton to fermions exists, provided that the relevant coupling constant 
is somewhat suppressed. If these fermions are identified with right-handed neutrinos, baryogenesis 
via non-thermal leptogenesis 1391 is, in principle, possible - in the case of HI, baryogenesis can be 
also accomplished if only the waterfall field is coupled to right-handed neutrinos. Note that, in our 
framework, the decay of the inflaton to right-handed neutrinos is also possible due to curvature-induced 
11241 couplings. However, the resulting decay width is reduced [24] compared to this given by Eq. ( 12.51) 
and so, the produced lepton asymmetry is lower than the expectations for all possible masses of right- 
handed neutrinos. On the other hand, since baryogenesis can be realized in a variety of ways - see, 
e.g., Ref. |[25ll40l - we opted not to complicate our presentation with secondary mechanisms which 
may or may not affect the inflationary observable quantities. 

It would be interesting to investigate if a similar realization of non-MI can be accomplished in the 
framework of SUGRA, along the lines of Ref. BP . In such a case, the inflaton of non-MCI could be 
identified with one of the right-handed sneutrinos. On the other hand, a possible SUSY version II15II16I1 
of non-MHI could become compatible with larger (and more natural) values of the relevant coupling 
constant k = A. 
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